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Matching of the Heavy-Light Currents with NRQCD Heavy and 
Improved Naive Light Quarks 

E.Gulez a , J.Shigemitsu a , M.Wingate b 

a Physics Department, The Ohio State University, Columbus, OH 43210, USA. 
institute for Nuclear Theory, University of Washington, Seattle, WA 98195, USA. 

One-loop matching of heavy-light currents is carried out for a highly improved lattice action, including the 
effects of mixings with dimension 4 0(1/ M) and O(a) operators. We use the NRQCD action for heavy quarks, 
the Asqtad improved naive action for light quarks, and the Symanzik improved glue action. These results are being 
used in recent heavy meson decay constant and semileptonic form factor calculations on the MILC dynamical 
configurations. 



1. Introduction 

Matching of the heavy-light currents between 
lattice simulations and full continuum theory is 
important for studies of heavy meson leptonic 
and semileptonic decays. As statistical and dis- 
cretization errors are reduced, chiral extrapo- 
lation and operator matching errors dominate, 
and more accurate matching calculations become 
necessary. As part of the work of studying 
heavy mesons, one-loop perturbative matching 
of the heavy-light currents was performed. In 
this talk we focus on the matching of the tem- 
poral currents Vo and A a . Our calculations are 
correct through 0(a s ,aa s ,a s /aM, ck s Aqcd/M)- 
Matching of the spatial components Vk and 
to the same level of accuracy has now also been 
completed. 
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2. Actions Used and Calculational Strate- 
gies 

The actions used are highly improved lattice 
actions. The AsqTad improved naive light quark 
action, which is correct through 0(a 2 ) was used. 
For the heavy quarks, we used an improved 
NRQCD action. This action is correct through 
(D(l/M 2 , a 2 ), and is the same action that is being 
used in numerical simulations on MILC dynami- 
cal configurations. For the gluon, the Symanzik 
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Figure 1. Comparison of the kinetic and per- 
turbative meson masses. The rightmost point at 
aMo — 2.8 corresponds to the physical B s meson. 
The data for this plot is from pQ . The points rep- 
resenting Mkin were shifted to the left for clarity. 
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Figure 2. Some results for the perturbative 
matching coefficients pj . The data for this plot is 
from 0. 



Figure 3. The effects of matching on fn s \/ M h s ■ 
At the value of the inverse B s mass, the effect of 
the perturbative correction is very small (within 
3%). The data for this plot is in pQ. 



improved glue action was used. 

Two independent methods were employed for 
the matching calculation:First, Mathematica was 
used to do the Dirac algebra and where neces- 
sary, to take derivatives with respect to exter- 
nal momenta. The result was converted into For- 
tran format and numerical integration was done 
with the Fortran Vegas program. Second, a C++ 
code originally developed for a similar calculation 
with other types of actions by C. Morningstar was 
modified and used. The derivatives were taken by 
automatic differentiation. 

3. Self Energy 

For massless Asqtad Quarks: 

Z q = l + a^C 1 / + C q ] + 0{a 2 s ) 
C\ R = J- [! + (£-!)] ln(a 2 A 2 ) 

O 7T 

q — °g + ° ? ' q 



u = +0(a 2 s ). 

(2) 

Mq = 0.767 for the plaquette definition of uq 

(2) 

and Uq = 0.750 for the Landau link definition, 
and A is the gluon mass. 
For heavy Quarks: 

Zq = 1 + a s [C Q R + C Q ] + 0{a 2 s ) 
Z M = 1 + a s C M + 0(a 2 s ) 
aE = a s CEo + 0(a 2 s ) 

C Q R = ^-[-2+(C-l)] Ma 2 A 2 ) 

Z q was previously calculated by H. Trottier using 
twisted boundary conditions as IR regulator. Our 
results for the lattice-continuum difference are in 
agreement. 
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4. Kinetic and Perturbative Masses for 
Heavy-Light Mesons 

Kinetic mass is obtained nonperturbatively 
from finite momentum correlators using the me- 
son dispersion relation; 
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where A = E(p) — E(0). The perturbative mass 
is estimated using the zero momentum correlator 
and the one-loop heavy quark mass renormaliza- 
tion 
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As shown in Fig. f , the results are in agreement 
with each other. 

5. Current Matching 

The lattice currents that are matched with the 
corresponding continuum ones are listed in 
The continuum currents are a linear combination 
of lattice ones: 



(A^CONT = Y, C t^ J A»)LAT, 



CONT 



6. Summary 

We have performed one-loop matching of the 
lattice NRQCD+AsqTad heavy-light current to 
the corresponding continuum current for different 
heavy quark masses. The calculation is correct 
through G(a s /aM). The perturbative matching 
coefficients are found to be small and well be- 
haved. These matching coefficients are important 
for the study of leptonic and semilcptonic decays 
of B and D mesons and of neutral B mixing. The 
results of this work is used for calculating the de- 
cay constants for B s , D s pQ and B These re- 
sults are also important ingredients in calculation 
of form factors /o, /+ 0. Similar perturbative 
matching for four fermion operators is underway. 
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The matching coefficients for the axial vector and 
vector coefficients are the same due to the good 
chiral properties of the AsqTad action; Cj fl = 

C. ** . The perturbative expansion of the matching 
coefficients goes as, 



Co,i — 
C 2 = 



1 + 
ap 2 



ap ,i + 0{a 2 
+ 0{a 2 ), 



The goal is to find the perturbative matching co- 
efficients pj. As seen in Fig. 2, these match- 
ing coefficients are small and well behaved, which 
shows that perturbation theory is working well. 
For the numerical values of the coefficients, see 
[3] .As an example, the effect of the matching on 
$ff s = Sh, \fM Ha is plotted in Fig 3. The differ- 
ence is within 3% around the physical B s . 



